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A new approach to investigate dineutron correlation 
and its application to ^°Be 

Fumiharu KOBAYASHi and Yoshiko Kanada-En'yo 
Department of Physics, Kyoto University, Kyoto 606-8502, Japan 

We propose a new framework by means of the dineutron condensate (DC) wave function 
to describe the dineutron correlation, which is characterized by the spatially strong corre- 
lation of a spin-zero neutron-neutron pair, in neutron-rich nuclei with an active deformed 
core surrounded by valence neutrons. Using the DC wave function for a 2a+2n system, 
which corresponds to a toy model for the ^°Be system, we investigate the neutron-neutron 
, correlation around the 2a core and discuss the mechanism of the dineutron formation at 

■ the surface of finite nuclei. To investigate dineutron correlations in realistic nuclear systems, 

we superpose the antisymmetrized molecular dynamics (AMD) wave functions and the DC 
I wave functions. Applying the AMD-I-DC method to ^"^Be, we show effects of the DC wave 

. functions in the ground and excited 0"'' states of ^"Be and discuss the dineutron correlation 

■^1^ ' in them. 

, §1. Introduction 

Recently, neutron-rich unstable nuclei are investigated intensively and various 
exotic phenomena have been discovered. The dineutron correlation is one of the 
Q I key issues attracting a great interest in physics of unstable nuclei. The dineutron 

^ ■ correlation means the strong spatial correlation between two neutrons in the spin- 

^ , , zero channel. Even though two neutrons do not form a bound state in free space, 

they attract each other rather strongly. It was suggested that the spatial correlations 
^ between two neutrons enhance in a low-density infinite nuclear matter,'^^^ in the 

, neutron-skin structure of medium-heavy nucleP^'^SJ and also in the halo structure 

On I of light nuclei.'^'"'^ In such nuclear systems, the strongly correlated neutron-pair 

' with a rather compact size can be regarded as a semi-bound or a virtually bound 

• , state called a dineutron. The dineutron correlation is considered to be an important 

^ I aspect to understand the structures of unstable nuclei. 

In studies of the asymmetric nuclear matter, it was shown that the strength of 
the two-neutron correlation changes as a function of the matter density .'^^^''^"^ 
Matsuo suggested the continuous transitions from a weak BCS-like coupling to a 
strong boson-like coupling in accordance with a decrease of the matter densitj^ (the 
BCS-BEC crossover, such as in a superfiuid Fermi gas^). The bosonic behavior of 
dineutrons is often discussed in relation to the a condensation suggested in a low- 
density symmetric nuclear matter.^® 

The dineutron correlation in neutron-rich nuclei has been investigated well so far. 
In the medium-heavy nuclei, Matsuo et al. have investigated the neutron-neutron 
correlation in O, Ca, Ni isotopes.® They used the HFB theory to deal with the 
correlation between two neutrons in those nuclei, and concluded that two neutrons 
in the neutron-skin region show the spatial correlation. Also in light-mass nuclei 
such as ^He, the dineutron correlation has been theoretically studied.'^'^ 
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For the neutron-halo nuclei, such as He and Li, where two valence neutrons 
are weakly bound to a core nucleus, the dineutron component was explicitly in- 
vestigated by means of three-body models assuming a core and two valence neu- 
tronsP^Sl''^^^ In Ref. [Ti|) . Hagino et al. suggested that the pair of neutrons 
having the strong angular correlation is distributed in the halo region of ^He and 
^^Li. In addition, they also have discussed the change in the correlation between 
two valence neutrons as a function of the distance from the core, associating with 
the BCS-BEC crossover in finite nuclear systems. ^"^^ Since most of the three-body 
models, however, have such an assumption as a spherical and inert core, they are not 
able to deal with systems having a deformed core and their applications are limited 
to some specific nuclei. 

In order to elucidate the properties and mechanisms of the dineutrons and their 
condensate states, it is significant to investigate dineutron behaviors in neutron-rich 
nuclei systematically. We, accordingly, have constructed a new model which can be 
applied to various nuclei having deformed structures to draw the universal properties 
of dineutrons at the nuclear surface. 

In this model, the nuclear system is assumed to consist of a deformed core 
and surrounding dineutrons. The deformed core is described with antisymmetrized 
molecular dynamics (AMD),'2SIi^23 and several pairs of neutrons around the core are 
written by the condensate wave function, in which dineutrons are moving in the S- 
orbit as as in the a condensate wave function.-^- We call this model "the dineutron 
condensate (DC) wave function". The AMD is a useful model to describe structures 
of light nuclei, such as cluster structures and deformations. By using the AMD to 
describe the core structure, the DC wave function, therefore, can contain the cluster 
structure in the core nucleus as well as the core deformation. Although the DC wave 
function is expected to be useful to analyze dineutron correlations in general nuclei, 
it may be too simple to describe details of realistic nuclear states because of the 
assumption of explicit dineutrons. Thus, besides the DC wave function, we adopt the 
AMD wave function for the total system, and superpose the AMD wave function and 
the DC wave function to construct "the AMD+DC wave function" . This AMD+DC 
method is a beyond AMD approach, which can efficiently supplement the AMD 
wave functions with dineutron correlation components. Since the AMD+DC method 
can be applied to general nuclei without specific assumptions, it may be useful to 
investigate systematically dineutron correlations in various nuclei. 

As a first step, we formulate the DC wave function and apply this framework to 
the system containing a 2a core and one dineutron, which corresponds to a toy model 
for the ^'^Be system. By calculating the energy expectation value of the dineutron 
with the DC wave function and analyzing its dependence on the dineutron size, 
we investigate the properties of a dineutron in the 2a+2n system and clarify the 
mechanism of the dineutron correlation enhancement at the nuclear surface. Then, 
superposing the 2a+2n DC wave functions and the ^''Be AMD wave functions to 
construct the ^°Be AMD+DC wave functions, we analyze the contribution of the 
dineutron components in 0"*" states of ^'^Be. Effects of the core excitation on the 
dineutron correlation in ^'^Be(O^) are also discussed. 

This paper is organaized as follows. In §2, we describe the framework of the 
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new approach. We explain the framework of the AMD wave function, the DC wave 
function, and the AMD+DC method, and describe the present method for the ap- 
pHcation to ^'^Be. In §3, we refer to the effective Hamiltonian used in the present 
calculation. In §4, we analyze the results of the 2a+2n DC wave function, and in 
§5, we show the results of ^"^Be calculated with the AMD+DC method and discuss 
the dineutron components in ^''Be. In §6, a summary and an outlook are given. 

§2. Framework 

We use two types of wave function to describe nuclear systems; the AMD wave 
function and the DC wave function. At first, we explain the formalism of each wave 
function. Then, we describe the present DC wave function for the 2a+2n system, 
and explain the AMD+DC method for the ^^Be system. 

2.1. AMD wave function 

The AMD method has an advantage that the nuclear structure can be repre- 
sented without specific assumptions such as axial symmetry and existence of clusters. 
In the AMD, the wave function |<?amd) for a nuclear system with the mass number 
A is described by a Slater determinant of A single-particle wave functions as 

|^AMd(-^)) = ^det[|^,)] (Z = {ZJ = {y„|J (* = 1,... 



= __det[|Vy,)|x^,)|r.)], (1) 

spatial part : (rlipY,) = (^-^^ exp[-i^i(r - Yif], (2) 

spin part : \Xi^) = CitI t) + i)' (3) 

isospin part : |rj) = \p) or |n), (4) 



where the i-th single-particle wave function is composed of the spatial part 
{ipYi), the spin part Ix^^) and the isospin part |rj). The subscripts of the spatial and 
spin parts, represent the variational parameters for each single-particle wave 

function. The spatial part is approximated by a Gaussian wave packet which spreads 
around the center Yi with the width l/^/2l^. The present AMD wave function is 
similar to that used in the fermionic molecular dynamics (FMD).'33'I^ In principle, 
the width parameters can be different for each Gaussian wave packet as Ref. |3^ . 
However, in the application to ^''Be in the present work, we use the common width 
u = 0.235 fm~^ for all single-particle wave functions for simplicity. 

2.2. DC wave function 

We here mention the general form of the DC wave function to deal with the 
dineutron condensate state consisting of a core and surrounding dineutrons. Later, 
we consider the specific case of one dineutron around the 2a core. 

In a DC wave function, a nucleus consists of a deformed core described with an 
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AMD wave function and some dineutrons around the core in the lowest S-orbit to 
form the condensate state. We formulate the general form of the DC wave function 
I^Dc)) describing an ^-nucleon system including fco-dineutrons and the core with 
A — 2kQ nucleons as follows. 



\<Pbc) = ^ n J d^Ynkexp 



Y2 y2 y2 

^ nkx _ nky _ ^ nkz 
52 52 52 



k=l ■ 

X det [l^/ii) • • • |'0A-2fco)l^nlt>IV'nl4,) " " " \i^nko^)\i^nko\)] , (5) 

where the single-particle wave functions of the nucleons composing the core, {i = 
1, • • • ,A — 2ko), are written by Gaussian wave packets in the same form as those in 
the AMD wave function described in Eq. ([1]) . They are parametrized by the Gaussian 
centers and the spin orientations, {Yi,^^} {i = 1, - ■ ■ ,A — 2ko). The wave function 
for the k-th. pair of neutrons, \ipnkt)\''Pnkl) {k = " > ^o) in the Slater determinant 
represents two neutrons with spin f and I, whose spatial parts are expressed by 
Gaussian wave packets as 

\lpnkx) = \'^nk)\x)\n), (6) 

spatial part : (rjV'nfc) = (--^^ exj)[-Unk{r - Ynkf], (7) 

1 

I'nk 



Kk 



spin part : |x) = I t>or| \). (8) 

Two neutrons have the common spatial wave functions with the width parameter 
Vnk and the Gaussian center Ynk-, that is, they spread around the point Ynk by an 
expansion h^k = l/\/2i/„fc. We use the parameter hnk instead of Unk hereafter. By 
integrating the Slater determinants in the weight of the Gaussian distribution with 
respect to the centers of dineutrons Ynk^ each dineutron spreads around the core 
in the S-orbit. We can calculate analytically the Gaussian integral for the energy 
expectation value of the Hamiltonian given in §3, and also the expectation values 
of the root mean square radius and density shown in §5. The parameter Bna in 
Eq.Q corresponds to the width of the dineutron distribution, and it is taken to 
be a common value for all dineutrons to describe a dineutron condensate state. In 
principle, Bna can be different values with respect to a = x, y, z. We, however, use 
the same value Bn = Bnx = Bny = Bm in the present calculations for simplicity. 
Then, by definition, the DC wave function contains dineutrons which correspond 
to the spin-zero neutron-neutron pairs with the relative orbital angular momentum 
I = and moving in the S'-wave around the core. The general DC wave function 
defined above has an analogy to the a condensate wave function proposed in Ref. [3T]) 
except for the existence of the core surrounded by dineutrons. As stated later, we 
need to project the DC wave functions to the angular momentum eigenstates for the 
description of realistic states. Note that the angular momentum projection of the 
DC wave functions is equivalent to that of the core wave functions, so the angular 
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momentum projection can be achieved by superposing the core wave functions. We 
also should stress that the DC wave function is a microscopic ^-body wave function 
with the full antisymmetrization between all nucleons including valence neutrons and 
core nucleons. 

2.3. DC wave function for 2a+2n 

Let us consider the particular case, ^''Be (2a+2n), where the system consists 
of two as and one dineutron (not a dineutron condensate state). As indicated in 
previous studies, for instance Ref. [35j) . most states of ^''Be have a rather developed 
2q core, so we assume the core as 2a. Below we explain the DC wave function for 
the 2a+2n system in detail. The 2Q+2n DC wave function is written as 



1^, 



'2a+2n/ 



1 



d Yn exp 



■ E 

cr=x,y,z 



y2 

^ no- 



det[[|^„l)]4[|V;a2)]1V;nt)IV;n4>]: 



(9) 



where \ipa) represents four wave functions composing an a, that is, pf, pi, nf and n|. 
We can write the ten-body wave function in the Slater determinant with a product 
of the spatial and spin-isospin parts, ^2a+2n X2a+2n, and then we rewrite the 

|^^2a+2n) aS folloWS. 



1 



2«+2n/ 



00! 



52 



X2a+2n 



(10) 



2a+2n 



{d, Yn, bn) = [i'al{ri) ■ ■ ■ ^^01(^4) ^a2(r'5) " " " ^02(^^8) V'n (^"9) V^n (^-lo)] 



Here we fixed the distance between two as to d and locate one a at (0, 0, d/2) and the 
other at (0,0,— d/2). For our purpose, we integrate Yn with the Gaussian weight 
and rewrite two valence neutron wave functions in Eg. pip with the relative and 
center of mass wave functions of two neutrons in a dineutron, 'ip2n and ipc- Then, 
the spatial part of Eq. ()10p can be rewritten as 



d Yn exp 
1 



y2 



'^2a+2nid, Yn, bn) 



[V'al('^l) • • • i^alir^) i'a2{r5) ■ ■ ■ V'«2(t'8) ^^271 (^')V'G (^g)] , 





^2 


(^) 





ipG = exp 



/32 



rg - rio, 
rg + rio 



(12) 
(13) 



P^^Bl + bl (14) 



The coordinates r and are the relative and center of mass coordinates for two 
neutrons. As seen in Eq. p3p and Eq.p4p. the relative and center of mass wave 
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functions are also written as the Gaussian forms whose width parameters are 6„ 
and /3, respectively, which satisfy the relation 6^ < /3 by definition. It means that 
two neutrons spread by the expansion b^i each other and the center of mass of the 
dineutron is distributed around the core by the extension /3. So we consider the 
parameter 6„ as the size of the dineutron, and f3 as the spread of the dineutron from 
the core, so that we characterize the behavior of the dineutron by 6„, and f3 instead 
of Bn- As a result, the 2a+2n DC wave function is specified by the parameters d, 
bn, and /3. 

The wave function V'2nV'G indicates the two- neutron wave function in the spatial 
part of the DC wave function before the antisymmetrization. In the case of bn ~ /3, 
the single-neutron wave function for two neutrons becomes a Gaussian located 
at the center of the total system (Bn ~ 0), and the two-neutron state corresponds to 
the uncorrelated limit of two valence neutrons moving around the core in the s-orbit 
with the width 5„ ~ /3. In contrast, when bn is relatively small compared with /3, 
two neutrons form a compact dineutron moving in the S-wave around the core. 

2.4. AMD+DC wave function for ^^Be 

In the DC wave function, all valence neutrons are assumed to couple to be spin- 
zero dineutrons. This assumption is not necessarily appropriate for realistic nuclear 
systems. In the present study, to investigate dineutron components in ^'^Be, we con- 
struct the AMD+DC wave function by superposing the AMD wave function and the 
DC wave function defined above. We use the AMD wave functions for ^'^Be(0]'') and 
^'^Be(O^) obtained by the method of variation after parity and angular momentum 
projections (VAP). The details of the VAP method in the AMD framework is de- 
scribed, for example, in Refs. [35]) .[36 ]) . We consider the superposition of the states 
|^AMD,fc) given by the AMD wave functions ^p^md{Z^^^), and the states |^DC,fc) 
given by the DC wave functions having various sets of the parameters (dk, (3k,bnk)', 
I^DC.fc) = \^i)c{dk, Pki bnk))- The numbers of the basis AMD wave functions and the 
DC wave functions are /camd and /cdCj a-^d k is the label for the AMD and DC wave 
functions, i.e., /c = 1, • • • , /camd for |<?AMD,fc) and k = /camd + !,••• , ^amd + knc 
for I^DC.fc)- We project those /camd + ^DC basis wave functions to the parity and 
angular momentum eigenstates with the projection operator 'P/^ which projects to 
the eigenstate of the total angular momentum J, the z-component of the angular 
momentum in the laboratory frame M and that of the body- fixed frame K, and the 
parity it. For 0^ states, the AMD+DC wave function is given by the superposition 
of the projected states as follows. 

^AMD *;amd+*;dc 

|!^AMD+DC> = ^ /fcT^Oo'l^AMD.fe) + Yl /fc^00^I^Dc(4, A, ^nfc)), (15) 

k=l fc=fcAMD + l 

The coefficients fk are determined by diagonalization of the Hamiltonian and norm 
matrices by solving the Hill- Wheeler equation,l33 



(16) 
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We calculate the angular momentum projection numerically by replacing the inte- 
grals of Euler angles with the sum of minute meshes. 

§3. Effective Hamiltonian 



The Hamiltonian used in the present calculation consists of the kinetic term and 
the effective nuclear forces, 



H = T-Tg + Fcent + Vl^S + Vconh 



(17) 



where T and Tq are the total and center of mass kinetic energy. In our framework, 
it is difficult to exactly remove the center of mass motion from the wave function 
because different width parameters are used in the DC wave functions. We, therefore, 
approximately treat the center of mass effect in the energy by subtracting Tq from 
T in the Hamiltonian. The potential term T4ent; ^LS and Vcoui are the central, 
spin-orbit and Coulomb force. We use the Volkov No. 2 interaction as V^ent!® 



Vr 



cent 




X X. 



HVr-MVaVr, 



(18) 



(19) 



where the operators Va and Vr are the spin- and isospin-interchanging operators. 
We fixed the parameters as vi = —60.65 MeV, V2 = 61.64 MeV, ai = 1.80 fm, 02 = 
1.01 fm. As for the spin-orbit force, we use the G3RS interaction,'22J'S3 



Vk exp 



k=l 



i<j 



a-k 



Vi^O) Lij ■ S 



(20) 



where Lij and Sij are the relative orbital angular momentum and the total spin of 
two particles, and the parameters are ai = 0.447 fm, a2 = 0.6 fm. The operator 
V(^0) projects two particles to the triplet-odd state. The Coulomb force Vcoui is 
approximated by a sum of seven Gaussians. 

In the present calculation, we adopt the same interaction parameters as those 
used in Ref. iU, that is, B = H = 0.125, M = 0.60 in Eq.^, and vi = -V2 = 
1600 MeV in Eq.dSO]). 



§4. A dineutron around a 2a core 

As the first step, we apply the DC wave function to the 2a-|-2n system to under- 
stand the properties of a dineutron around the 2a core. We discuss the mechanism 
of the dineutron correlation enhancement by analyzing the energy of the dineutron 
as a function of the dineutron size bn and the spread of the dineutron distribution 
/3. 

We define the dineutron energy, , measured from the limit where the 2a core 
exists at the origin and two neutrons are free, i.e. in the zero-momentum state. 



E2nid,(3,bn) = E2a+2nid,(3,bn) - E2a{d), 



(21) 
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2ct+2n {d,/3,bn)) 

^2a+2n{a, P,On = —TT / , o , / , a , x\ ; {^^) 

(^2a+2n(a, P, On ) |^2a+2n (a, P, bn)) 

(23) 

where E2a{d) is the energy expectation value for the DC wave function in the infinite 
bn and /3 hmit. With a fixed d for the 2a distance and a fixed /? for the expansion 
from the core, the energy E2nid, /3,bn) is regarded as the binding energy of two 
neutrons virtually confined around the 2a core. 

At first, we consider the behavior of the dineutron energy E2n when its size and 
expansion from the core, bn and /3, are changed. Here we fixed the a-a distance d 
to 1, 2, 3, 4 fm. We show in Fig. [Jthe dineutron energy E2n plotted as functions of 
{d, /3, bn)- Since the qualitative feature does not depend on the parameter d so much, 
we focus only on the case d = 2 f m (the upper-right figure in Fig. [T|), which gives 
the lowest solution for the energy E2a+2n of the total 2a+2n system. The energy 
minimum exists at {f3,bn) ~ (1.9,1.5). As the viewpoint of the dependence on the 
dineutron size bn, the energy increases rapidly when bn deviates from the minimum 
point bn ~ 1-5 fm. This means that a dineutron around the 2a core favors a fairly 
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B (fm) 

3 4 
b„ (fm) 

Fig. 2. The dineutron energy, E2n, in the 2a+2n system for the fixed 2a distance d = 2 fm. The 
left figure is the energy surface on the /3-6„ plane same as the figure of d = 2 fm in Fig. [T] The 
right one is i?2n plotted as a function of bn for the fixed /3 = 3, 5, 7 fm, which corresponds to 
the section on the line /3 = 3, 5, 7 fm of the energy surface shown in the left. 



compact size. On the other hand, the energy increases only gradually when the 
spread from the core /? increases from /? ~ 1.9 fm. Thus, the dineutron can extend 
far from the 2a core to some extent. These results suggest that the dineutron around 
the 2a core tends to be distributed widely keeping a compact size. 

Subsequently, we discuss the origin of the compact dineutron formation below. 
We analyze the 6„-dependence of the dineutron energy with a fixed /3-value in more 
detail. The DC wave function with a fixed f3 corresponds to the state where a size- 
chageable dineutron is confined to a finite region around the core. In Fig. [21 we 
plot the dineutron energy as a function of the dineutron size 5„ with /3 fixed to 
be 3, 5, 7 fm. When /3 = 3 fm, that is, the dineturon is confined relatively close 
to the core, in other words, near the surface of the nucleus, the energy structure 
shows a high barrier at the point 6„ ~ /3 and a deep pocket at bn ~ 1.5 fm. Such 
structure disappears gradually as the spread from the core, /3, increases. The state 
6„ ~ /3 corresponds to the uncorrelated limit of two valence neutrons, while the state 
bn ~ 1.5 fm corresponds to the state where two neutrons correlate strongly to some 
extent to form a compact dineutron. When two valence neutrons are distributed 
with /3 = 3 fm, i.e. near the nuclear surface, they are strongly correlated to be a 
rather compact size due to the pocket in their energy, and the high barrier prevents 
two neutrons from separating each other by a large distance. 

As mentioned above, the existence of the barrier is essential to enhance the 
dineutron correlation. This barrier structure suggests that two neutrons in the un- 
correlated limit feel the effective repulsion. In order to clarify this mechanism of 
developing the dineutron, we decompose the energy E2n into the kinetic part and 
the potential part (Fig. [3]). The potential energy decreases monotonously as (3 or 
bn increases, as can be interpreted easily. As the dineutron size bn becomes small, 
two neutrons gain the potential energy due to the attraction of the neutron-neutron 
interaction. The decrease of /3 corresponds to the shrinkage of the confining region 
of the dineutron, and therefore the attraction from the core increases. This simple 
6,1-dependence of the potential term does not contribute to the barrier structure of 
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the dineutron energy i?2n- 

On the other hand, the kinetic energy has a distinctive structure. The kinetic 
part of the dineutron energy for /3 = 3 fm, where the dineutron are distributed 
near the nuclear surface, shows the plateau at the region ~ /3, while that for 
/3 = 5, 7 fm decreases monotonously as bn increases. This monotonous decrease 
is due to the uncertainty principle for two neutrons in the dineutron. As the size 
of a dineutron bn becomes smaller, the uncertainty principle enlarges the kinetic 
energy of two neutrons in it, and vice versa. When the dineutron is distributed 
near the surface of the nucleus, in addition to such an effect, another effect arises 
from the Pauli blocking against the core, and it is the main origin of the plateau. 
When the dineutron is distributed near the core, the Pauli exclusion principle works 
significantly between the neutrons in the dineturon and those in the core. The state 
at 6n ~ /3 corresponds to the uncorrelated limit state where each valence neutron 
occupies the isotropic s-orbit with the width bn around the core. Since the neutrons 
in the core already occupy the s-orbit with the width l/\/2^', two valence neutrons 
feel the strong Pauli repulsion from the core neutrons. The increase of the Pauli 
repulsion constructs the plateau at the region bn ^ /5 in the kinetic term. Whereas, 
when the dineutron size 6„ is small relatively to /?, two neutrons have the strong 
angular correlation and their single-particle wave functions contain components of 
higher angular momentum orbits. Then the Pauli blocking effect from the core 
becomes weak, and the Pauli repulsion in the kinetic term is suppressed. In this 
way, the plateau of the kinetic part of the dineutron energy for /? = 3 fm originates 
in the Pauli repulsion from the core. In cases of /? = 5 and 7 fm, since the dineutron 
distributes far from the core and the Pauli repulsion is generally weak, there is no 
clear structure of such a plateau in the kinetic energy. 

In the above analysis, it is concluded that the high barrier and the deep pocket 
seen in the dineutron energy E2n for /3 = 3 fm are constructed by the plateau 
structure of the kinetic part and the attraction of the potential part. The effects 
from the core are essential for these structures, i.e. the Pauli repulsion from the core 
yields the barrier, and the additional attraction from the core contributes to the 
deep pocket. Because of the barrier and pocket structure in the dineutron energy, a 
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dineutron at the nuclear surface favors a suitable compact size. Also in the a+2n 
and 3a+2n, we find a similar mechanism for dineutron formation, though the depth 
of the energy pocket depends on the core. Since the dineutron size at the energy 
pocket is common for these nuclei, it turns out that the dineutron correlation occurs 
in the circumstance of the s-orbit occupied by core neutrons and the unoccupied 
p-orbits. This mechanism of the dineutron correlation enhancement due to the Pauli 
repulsion from the core might occur in general nuclei. The favored size, however, 
may depend on the allowed orbits for each nucleus. 

§5. Contribution of a dineutron component to ^°Be states 

5.1. Superposition of AMD and DC wave functions 

In this section, we apply the AMD+DC method to ^^Be and analyze the dineu- 
tron component in ^^Be states. Since, here, each DC wave function has an isotropic 
parameter Bn (= Bnx = Bny = Bnz), we discuss only 0^ states which are thought 
to be described properly in this parameter fixing. We adopt the AMD wave func- 
tions corresponding to the ground state and the first 0"*" excited state obtained by 
the VAP. We superpose DC wave functions having various bn- and /3-values, with 
two AMD wave functions of ^'^Be(OJ['", 0^). We superpose DC wave functions not 
only with d = 2 fm, but also having different d-values. The superposition of DC 
wave functions with different d-values corresponds to the situation where two as 
can move relatively and, consequently, the consideration of the core excitation. We 
also perform the fixed-d calculations and compare the results with the superposing-d 
calculations to see effects of the core excitation. Thus, we consider the following 
superpositions, 



(i) !P' = ^/fcPOo+^AMD,A:, (24) 

fc=l 
18 

(ii) ^ = Y, fkV^o'^Dcid = 2, (3k, bnk), (25) 

fc=3 

2 18 

(iii) ^ = Y, An°o^^AMD,fc + fk'Poo'^^cid = 2, (3k, bnk), (26) 

k=l k=3 
2 66 

(iv) ^ = Y, fkr^0'^AMD,k + fk'P^0'^Dc{dk,^k, bnk). (27) 



k=l fc=3 



For the parameters of basis DC wave functions, we adopt /? = 2, 3, 4, 5 fm, and four 
bn-vslues for each (3 in the following way^^ 



(i-l)/(imax-l) 



bni = bnl X i^-^j (i = 1, • • • ^max), (28) 

where we choose imax = 4, bni = 1.5 fm and bni,^^^ = /5 — 0.2 fm. The parameter 
d is chosen to 1, 2, 3, 4 fm, when superposing different-d DC wave functions ((iv)). 
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so that the number of the basis DC wave functions in the superposition is A;dc = 16 
for (ii), (iii), and /cdc = 64 for (iv). We hereafter label the wave functions (i) AMD, 
(h) DC, (iii) AMD+DC((i = 2) and (iv) AMD+DC(Z'd). 

5.2. Energy of 0^ states in ^^Be 

First of all, we show the energy spectra of AMD, DC, AMD+DC((i = 2) and 
AMD+DC(i7rf) for comparison in Fig. [H The lowest and second lowest levels in each 
calculation, (i), (iii) and (iv), correspond to the 0^ and 0^ states and their ener- 
gies are -58.64, -50.03 (AMD), -59.27, -50.34 (AMD+DC(d = 2)), and -60.35, 
-50.93 (AMD+DC(Z'd)) in MeV, respectively. In addition to the 0^ and 0^ states, 
many states are obtained by the diagonalization of the AMD+DC in the higher en- 
ergy region. Although most of them show features of unbound continuum states, we 
obtain an indication of the possible candidate for a developed dineutron state in the 
fifth state in AMD+DC(i7^). We will refer to that state in the end of this section. 

In the AMD wave functions describing the 0^ and 0^, the energy of the 0^ is 
higher than that of the 0^ by 8.61 MeV (6.18 MeV for an experimental value). In 
studies of ^''Be such as Refs. [35ll and H3]) . it is indicated that the 0^ state has the 
2a surrounded by two neutrons distributed predominantly in the tt^ configuration 
in terms of the molecular orbit, while the 0^ state has a remarkably developed 2a 



-40 r 



-45 - 



-50 - 



-55 



-60 - 



-65 



AMD 



DC 



AMD-HDC(d=2) AMD-^DC(I:d) 



Fig. 4. The energy spectra of 0+ states of "Be. The labels AMD, DC, AMD+DC(d = 2) and 
AMD-|-DC(£'d) indicate the spectrum calculated with only AMD wave functions (i), the one 
with only DC wave functions fixed d to 2 fm (ii), the one with the AMD-I-DC wave function 
whose parameter d is fixed to 2 fm (iii) and the AMD-f DC wave function whose parameter d 
is superposed (iv), respectively. The calculated threshold to a+a+n+n (—55.23 MeV) is also 
shown. 
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core structure with two neutrons in the configuration. 

Let us consider the energy gain in the Of and 0^ contributed by the DC 
wave functions. Comparing the energy between AMD, AMD+DC((i = 2) and 
AMD+DC(i7rf), we find that some energy is gained by mixing DC wave functions. 
The energy gain in AMD+DC(i7d) (~ 1 MeV) is larger than that in AMD+DC(d = 
2) in both states. It indicates that the core excitation significantly contributes to 
at least the energy of those states. Not only the energy but also the dineutron 
component is affected by the core excitation as shown later. These results indicate 
that, in mixing the 2a+2n DC wave functions, the variation in the a-a distance is 
important, and generally, the explicit core structure, such as the core deformation 
or polarization, may enhance the dineutron component. 

We compare our results with the energy spectra of other calculations in Refs. l4ip . 
H3|) as well as the experimental ones in Fig. [5l We show the spectra calculated with 
the molecular orbital (MO) model in Ref. H3]) and those with the /3-7 constraint 
AMD in Ref. 43). In both works, the same effective interactions were adopted. The 
only difference in the used interactions is the strength of the spin-orbit force, the 
parameters vi = —V2 = 1600 MeV were used in Ref. HTj) as well as the present 
calculation, while vi = —V2 = 2000 MeV were adopted in Ref. H3|) . Both works 
suggested the Oj|~ state in addition to the Of and 0^ states. The present results 
of the 0^ excitation energy and the position of the 2a+2n threshold energy are 



20 r 



15 h ^ 



10 - n+ 0^; 



0^2 Oj_ 



0*^ 0*^ 0*^ 0*^ O'^ 



EXP. AMD+DC{d=2) AMD+DC(E(j) ltagaki[43] Suhara[41] 



Fig. 5. The energy spectra of 0^ states of ^"Be on the basis of the ground state 0^ energy. The one 
labeled with Itagaki[43] is the calculated values with the molecular orbital model in Ref. 43), and 
Suhara[41] is the values with the /3-7 constrained AMD in Ref. 141). The broken lines correspond 
to the threshold to a+a+n+n for the experiment and each calculation. 
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Table I. The root mean square radii of matter, proton and neutron of ^"Be (0^,0 J) calculated with 
the AMD, AMD+DC(d = 2) and AMD+DC(X'd) wave functions. The labels are the same as 
that in Fig. |4l The experimental value of the matter radius of the 0^ state is also shown. 





AMD 


AMD+DC(d = 2) 


AMD+DC(rd) 


EXP. 




2.28 


2.24 


2.32 


2.30 ± 0.02 


' rms 


2.22 


2.22 


2.31 




r" 

' rms 


2.32 


2.34 


2.41 




O2 ^rms 


2.84 


2.84 


2.94 




' rms 


2.72 


2.68 


2.66 




' rms 


2.93 


2.93 


3.14 





consistent with those of Ref. i41p by Suhara et al. They underestimate the relative 
position of the 2a + 2n threshold energy to the 0^ and 0^ states compared with 
the experimental data and also with the work in Ref. B3ll . The main reason for 
the difference in the threshold position between the present and Itagaki et aVs 
calculations is the difference in the spin-orbit interaction parameter. When we use 
the same interaction parameters as those used in Ref. I43p . we obtain the better result 
of the threshold position. For the quantitative reproduction of the energy levels as 
well as the threshold energy, we need fine tuning of the effective interactions. 

5.3. Radii and density for proton and neutorn distributions 

We calculate the root mean square (r.m.s.) radii of the 0^ and 0^ states to 
discuss the expansion of the neutron ditribution. The r.m.s. matter radius r^j^s for 
a normalized wave function is defined as 

1 ^ 

{r,^s? = {r^) = -Y,mr}\^). (29) 

i=l 

The ones of proton and neutron are defined in the same manner. We also calculate 
the proton and neutron density. The proton (neutron) density, Pp(n), defined as a 
function of the distance from the core is defined as follows. 

Z(N) 

Pp{n)ir) = (1^1 Yl '^(^ - ^^)^n,P{n)\'^)- (30) 
1=1 

The calculated r.m.s. radii are listed in Table. [H and the density is shown in Figs. [6] 
and [71 

In the 0^ state, the radii of AMD+DC(Z'^) increase relative to those of AMD. 
Since not only the radius of neutron but also that of proton does become large, the 
radii increases mainly due to the a cluster development. Comparing the neutron 
density of AMD, AMD+DC((i = 2) and AMD+DC(Z'd), shown in Fig. El it is seen 
that the neutron density in the inner region of the nucleus (r < 2 fm) is not so 
different between these wave functions. As seen in the r > 4 region, the expansion 
of the neutron density can be described by superposing the DC wave functions. 
However, the difference between the proton and neutron density is not so remarkable 
as the neutron-halo or skin structure. 
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'Be(O^) 




01 2345678 
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^°Be(0*,) ^°Be(0*,) 




012345678 012345678 
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Fig. 6. Top; the neutron density calculated with AMD, AMD+DC(d = 2), AMD+DC(rd). 
Bottom- left; the density of proton and neutron calculated with AMD. Bottom-right; those 
calculated with AMD-f DC(X'd). 

Subsequently, we discuss the radii and the density of the 0^ state. It can be seen 
that the 0^ state has a rather extensive structure compared with the 0^ state. It 
is interesting that the radius of neutron increases remarkably though that of proton 
decreases shghtly. We can see such an expansion of a neutron, also in the neutron 
density of ^°Be(02^) shown in Fig. [71 The expansive structure of the valence neutrons 
described by the mixture of DC wave functions is well represented also in comparison 
of the proton and neutron density (the bottom-right panel of Fig. [7]). The neutrons 
are distributed rather extensively compared with protons, different from the 0^ state. 
However, as we discussed later, this extensive neutron structure is not necessarily 
due to the components of strongly correlated two neutrons distributed to the outer 
region but due to mainly those of uncorrelated two neutrons. 

5.4. Dineutron component in ^^Be(0^) and ^^Be(02) 

In the AMD-I-DC calculations, the 0^ and 0^ states contain the dominant 
AMD wave functions with the mixture of DC wave functions. In the 0'^ state of 
AMD-hDC(<i = 2), the overlap with AMD is 97%, and the contribution of DC 
wave functions is minor. On the other hand, the AMD component in the 0^ state 
of AMD-hDC(Z'd) decreases to 89%. In the 0^ state, the AMD component in 
AMD-hDC(d = 2) and AMD-FDC(Z'd) is 98% and 82%, respectively. Therefore, 
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'Be(0*2) 




012345678 
r(fm) 



^°Be(0*2) ^°Be(0*2) 




r (fm) r (fm) 

Fig. 7. Top; the one-neutron density in the state calculated with AMD, AMD+DC(d = 2), 
AMD+DC(£'d). Bottom-left; the density of proton and neutron in the 0^ state calculated with 
AMD. Bottom-right; those calculated with AMD-fDC(rd). 

the mixture of the dineutron component described with the DC wave functions is 
enhanced by superposing various d-values, i.e. considering the 2a relative motion. 
However the AMD component is still dominant and we suppose that the main prop- 
erties of the AMD-I-DC states may not be so much different from those of the AMD 
states. Despite its a little mixing amplitude, the mixture of the DC component can 
significantly contribute to the energy of the system and the dineutron component far 
from the core. In the following, we analyze the dineutron component in these states 
in detail. 

To investigate the dineutron component, we calculate the overlap of the AMD-I- 
DC{d = 2) or AMD+T)C{U(i) wave function with the DC wave function whose 
parameters 6„, /3 and d are fixed to certain values. The overlap of a normalized wave 
function jlZ') with a DC wave function \$Dc{d, I3,bn)) is written as 

N{d,^,bn) = m'pDc{d,(3,bn))\^ (31) 

This quantity is an indicator which shows how the dineutron component contributes 
to the state. We denote the overlap of AMB+BC{d = 2) and that of AMB+BCiUd) 
with ^]:,c{d, (3,bn) as Nd=2{d, (3,bn) and Nx;^{d, I3,bn), respectively. The calculated 
values of N(i=2id, /3,bn) for d = 2 fm and those of Nx;^{d, I3,bn) for d = 2,3, 4 fm 
are plotted on the /3-6„ plane in Fig. [51 The panels in the left row of Fig. [5] show 
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(a) Nd,2 (d=2) (e) N^.j (d=2) 




1234567 1234567 



P (fm) P (fm) 

Fig. 8. (a)-(d); the overlap of the 0^ state with DC wave function whose parameter d is fixed (the 
fixed value is represented in the parenthesis) on the surface. The label Nd=2, or A^^^, 

indicates the overlap of the 0^ state obtained by the AMD+DC(d = 2), or the AMD+DC(X'd) 
wave function. They are explained in the text, (e)-(h); the overlap of the in the same manner 
as the 0^. 
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the overlaps of the OJ" state, and those in the right row show the overlaps of the O2 
state. Since we deal with the parameter bn as the size of the dineutron and /3 as 
the spread of the dineutron distribution from the core, we consider the amplitude in 
the small 6„ and relatively large f3 region as the component of a developed compact 
dineutron. 

At first, we discuss the dineutron component in the 0^ state. As shown in 
Fig. [8|^b)-(d), AMD+DC(Z'(i) contains a large component of the DC wave function 
whose d is 2, 3 fm. The rather large component near (/3, bn) ~ (2.0, 1.5) (the pocket 
of the dineutron energy, see Fig. [1]) corresponds to the rather compact dineutron 
at the nuclear surface. Seeing these overlaps, there is the tail characterized by 
the relatively gradual decrease from the peak (/3, 6„) ~ (2.0, 1.5) toward the large /3 
region along bn ~ 1.5 fm. This component is the dineutron distributed away from the 
core to some extent keeping a compact size. Interestingly, the significant amplitude 
of the compact dineutron with 6„ ~ 1.5 fm in the large /3 region can be seen even 
in Njjjid = 4) for the spatially developed two a clusters. However, the maximum 
overlap with the DC wave function having d = 4 fm is half of those for d = 2,3 fm, 
which means that the component of the extremely developed two a clusters is not 
dominant in the state, unlike the 0^. 

To see how far the compact dineutron is distributed from the core, we show 
the overlap N{d, f3,bn) with the fixed dineutron size 6„, = 1.5 fm in Fig. [9j In 
the left figure of Fig. H the overlaps, Nd=2 and N^^, of AMD+DC(d = 2) and 
AMD+DC(i7rf) are enhanced compared to that of AMD in the region /? > 2 fm, 
where the dineutron is distributed away from the core. The remarkable increase of 
the overlap in this region indicates that the mixture of the DC wave functions plays 
a significant role in describing the developed dineutron. Next, we compare Nd=2 
and Njj^ to see the effect of core excitation. The maximum of each overlap exists at 
/3 ~ 2 fm, which corresponds to the vicinity of the nuclear surface, and there is no 
pronounced difference between them near this point. That means the variability of 




123456789 123456789 
P (fm) |3 (fm) 



Fig. 9. Left; the overlap of the OJ" states with the DC wave function where 6„ = 1.5 fm and d — 2 
fm. The label A'^amd, Nd=2, or A^^^ indicates the overlap of the 0]^ states obtained by the 
AMD, AMD+DC(d = 2), or AMD+DC(X'd) wave function. Right; the overlap of the 0+ states 
described with the AMD+DC(X'ti) wave function with the DC wave function having 6„ fixed to 
1.5 fm and d to 2, 3, 4 fm. 
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the a-a distance does not affect the development of the dineutron correlation from 
the inside to the surface of ^"^Be. The difference is seen in the region /3 > 2 fm, where 
the dineutron is distributed distantly from the core to some extent. Nx:^ is a bit 
larger than N(i=2 there, so that this increase suggests that the core excitation, such 
as the relative motion of clusters, may contribute significantly to the enhancement 
of the dineutron correlation. 

We show the overlap Nx;^{d, f3, hn) with the fixed dineutron size 6„ = 1.5 fm for 
DC wave functions having various d-values in Fig. [9j The dineutron component near 
the nuclear surface (/? ~ 2 fm) is greater in the case d = 2, 3 fm than d = 4 fm. The 
value of (3 giving the maximum of Nx:^ becomes larger as the parameter d of the 
overlapped DC wave function goes from 2 fm to 4 fm. That means the dineutron 
develops spatially and becomes distributed away from the core in accordance with 
the development of two a clusters. 

Subsequently, we consider the dineutron component in ^''Be(O^). As seen in 
Fig- [Stf)-(h), the 0^ state has the largest overlap with the DC wave function having 
d = 4 fm at ~ (1.5,1.5), which corresponds to not the compact dineutron 

near the nuclear surface but two neutrons distributed in the middle of the distant 
two a clusters. This behavior is consistent with the molecular orbital structure with 
the (T^ configuration, which has the nodal structure along the axis 2a located on. 
In any cases of d = 2, 3 and 4 fm, the overlaps have not too small amplitudes in 
the almost whole region, where the distance between two neutrons and their spread 
from the core are large. Namely, in the present calculation, two valence neutrons in 
the 0^ state are distributed extensively from the core with less correlation. 

5.5. Candidate for the resonance state Oj^ 

In the preceding studies, the Oj|' state having a developed dineutron around two 
developed a clusters were predicted ,'^''331' though there is no experimental evidence 
yet. In the present calculation, many states besides the 0^ and 0^ states are ob- 
tained in the AMD+DC wave functions. Although most of them are the insignificant 
continuum states, there would be a resonance-like state with a developed dineutron 
component. In Fig. [TOl we show the overlap, iVi;^, of the AMD+DC(i7rf) with the 
DC wave function for typical examples corresponding to continuum states, and that 
for the candidate of a resonance state. We choose the parameter d = 4 for the over- 
lapped DC wave function, for the resonance state is thought to contain the clusters 
bound loosely each other so that the overlap should be large when two as are distant. 

The overlap of the third lowest state in AMD+DC(i7£;) shown in the top-left 
panel of Fig. [10] is the typical one for a continuum state. This state has an enormous 
overlap with the DC whose both parameters 6„. and /3 are large. This tendency 
reflects that the three-body unbound feature of two valence neutrons and the core. 
Another example is shown in the top-right of Fig. [10] for the forth lowest one in 
AMD+T)C{U(i)- Two peaks at bn ~ 1-7 and 5 fm and a node at 6„ ~ 3 fm are seen 
in this figure. This also shows a feature of a continuum state. 

The state shown in the bottom panel of Fig. [10] is the overlap of the fifth lowest 
state, which we expect to be the candidate for a resonance state. This state have a 



20 



KOBAYASHI and Kanada-En'yo 




N. (d=4) 



0.7 




5 r 


0.6 




4.5 - 


0.5 




4 - 


0.4 
0.3 


E 
xf 


3.5 - 
3 - 
2.5 - 


0.2 




2 - 


0.1 




1.5 - 







1 L 




|3(fm) 



0.4 
0.35 
0.3 
0.25 
0.2 
0.15 
0.1 
0.05 


1 2 3 4 5 6 7 

|3(fm) 

Fig. 10. Top-left; the overlap of the third lowest state in the AMD-l-DC(X'd) with the DC wave 
function having d — 4. Top-right; that of the forth lowest state. Bottom; that of the fifth lowest 
state, which is expected to be a resonance state. Note that the scales of amplitude in the figure 
of the forth and fifth states are different from so far. 



peak in the moderate /3 and bn region {13, bn) ~ (3,2), which represents a relatively 
smaller distance between two neutrons and that between the dineutron and the core 
than those in the cases of continuum states mentioned above. The excited energy 
of this state is close to the ones calculated with different frameworks in preceding 
studied^'133 (see Fig. [5]). We, however, cannot yet conclude that this fifth state is 
a resonance with only the evidence we have shown here. Since this state may be 
combined with contiuum states, the matter radius and density do not converge. So 
we need to analyze it more closely in a different manner in future. 

§6. Summary 

We have proposed a new approach to investigate the dineutron correlation. To 
incorporate the dineutron component, we introduce "the dineutron condensate wave 
function", for short, the DC wave function. The wave function can describe a system 
composed of a deformative core and surrounding dineutrons which condensate in 
the lowest S-orbit. In the AMD-I-DC method for the investigation of the dineutron 
correlation in neutron-rich nuclei, DC wave functions are combined with AMD wave 
functions. 

We have applied the DC wave function to a particular system which consists of 
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a 2a core and one dineutron. Analyzing the dineutron size dependence of the energy 
of the 2a+2n system, we discuss the mechanism of the dineutron formation around 
the 2a core. We have found that the dineutron favors a relatively compact size when 
distributed near the nuclear surface because of the Pauli blocking effect from the 
core. Core effects are essential for the dineutron formation, and this mechanism may 
lead the dineutron correlation enhancement in general nuclei. 

We have investigated the dineutron correlation in O"*" states of ^''Be with the 
AMD+DC method. The obtained wave functions for the Of and 0^ states have the 
dominant AMD component with slightly mixing of the DC component. The mixing 
component of the DC wave functions enhances when we take into account the core 
excitation by superposing DC wave functions with different a-a distances. In the 
0^ state, we have found a tail structure of the dineutron in the region far from the 
core. In spite of the dineutron tail, the component of the strongly correlated two 
neutrons does not contribute so much to such one-body observables as the neutron 
radius and density, unlike halo nuclei. It may be because two neutrons are somehow 
strongly bound in the ^*^Be system. 

In addition to analyzing the dineutron component in the Of and 0^ states, we 
have searched for the 0^ state having a developed dineutron, which were predicted 
theoretically by preceding studies .SD'^S We have found a candidate for the state in 
the present results. We, however, could not conclude that the candidate state is a 
resonance because it contains the component of continuum states and the radius and 
density do not converge. We need further investigations by dealing with resonant 
and continuum states properly. 

In the present work, we have applied the new framework, "the DC wave func- 
tion", to a virtual system 2a-|-2n, and its practical application, "the AMD-I-DC 
method" , to a simple system ^'^Be. We still have a room for the improvement of this 
framework, for example adopting a deformed distribution of dineutrons. However, 
even in the simplest application of this work, we have obtained important conclusion 
for the mechanism of the dineutron formation and the possibility of the effect of the 
core structure to enhance the dineutron correlation. And also we have confirmed the 
utility of this framework to describe the dineutron correlation in the nucleus which 
contains a deformed and excited core. We expect that our framework can be applied 
to various nuclei without critical assumptions. 

In future, we will apply the AMD-I-DC method to various nuclei, in particular 
neutron-rich nuclei which are characterized with exotic structures such as the neutron 
halo and are expected to contain a developed dineutron. For instance, ^^Be is a 
challenging problem since it may be a two-neutron-halo nucleus with a deformed 
^^Be core. One of our final objectives is to clarify the universal properties of the 
dineutron correlation by systematic investigations of the dineutron correlation in 
diverse circumstances. Besides the case that the systems include one dineutron, 
we also would like to consider the system with a few dineutrons and analyze the 
property of a dineutron condensate state with developed dineutrons near the low- 
density surface. 
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